The origin of grand minima in the sunspot cycle 
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One of the most striking aspects of the 11-year sunspot cycle is that there have been times in 
the past when some cycles went missing, a most well-known example of this being the Maunder 
minimum during 1645 — 1715. Analyses of cosmogenic isotopes ( 14 C and 10 Be) indicated that there 
were about 27 grand minima in the last 11,000 yr, implying that about 2.7% of the solar cycles had 
conditions appropriate for forcing the Sun into grand minima. We address the question how grand 
minima are produced and specifically calculate the frequency of occurrence of grand minima from a 
theoretical dynamo model. We assume that fluctuations in the poloidal field generation mechanism 
and in the meridional circulation produce irregularities of sunspot cycles. Taking these fluctuations 
to be Gaussian and estimating the values of important parameters from the data of last 28 solar 
cycles, we show from our flux transport dynamo model that about 1-4% of the sunspot cycles may 
have conditions suitable for inducing grand minima. 

PACS numbers: 



A few years after the initiation of telescopic observa- 
tions of sunspots in 1610, there was a period from 1645 
to 1715 when very few sunspots appeared on the face of 
the Sun. This period is known as the Maunder minimum 
Although reliable sunspot data did not exist before 
1610, the solar activity at earlier times can be studied 
from the analyses of the abundances of cosmogenic iso- 
topes like 14 C in old tree rings and 10 Be in polar ice 
0| . When sunspots are absent, the magnetic field in the 
solar wind becomes weak and more galactic cosmic rays 
can reach the Earth, producing more of such radioactive 
isotopes in the atmosphere. Analyses of these isotopes 
indicate that there have been about 27 grand minima 
of different durations in the last 11,000 years Since 
there were about 1,000 solar cycles during this period, 
the occurrence of 27 grand minima implies that about 
2.7% cycles had conditions appropriate for forcing the 
Sun into grand minima. Furthermore this study showed 
that the Sun was in the grand minima state for about 
17% of the time. We also mention that there are evi- 
dences that some solar-like stars show grand minima Q . 
Therefore it is very important to understand the physics 
of the origin of the grand minima and the probability of 
occurrence of such grand minima. Moreover, sunspot cy- 
cles have an important effect on the space environment 
and the Earth climate system 0, Hj . Therefore, under- 
standing the grand minima is also important to the space 
weather and Earth climate communities. 

It appears that the most promising theoretical model 
for the solar cycle at t he p resent time is the flux trans- 
port dynamo model I9 H1 ll . which has been reviewed re- 
cently by Charbonneau [12| and Choudhuri fl3j | . Let us 
consider the question how irregularities arise in solar cy- 
cles. One important assumption of the flux transport 
dynamo model is that the poloidal field is produced by 
the Babcock-Leighton mechanism, in which tilted bipo- 
lar sunspots give rise to the poloidal field after their de- 



cay. The amount of poloidal field generated depends on 
the tilt angle of the bipolar sunspots. This tilt is pro- 
duced by the action of the Coriolis force acting on the 
magnetic flux tube rising through the solar convection 
zone due to magnetic buoyancy 14j and the average tilt 
at a solar latitude is given by Joy's law. However, there 
is a large scatter in the tilt angles around the average 
given by Joy's law 15-], presumably due to the effects 
of turbulence on the rising flux tubes (l6| . Hence the 
Babcock-Leighton mechanism has an inherent random- 
ness because of which the poloidal field generated at the 
end of a cycle would vary from one cycle to the other 
[l7j |. The second source of irregularities in solar cycles 
comes from the fluctuations in the meridional circula- 
tion of the Sun, which plays a crucial role in the flux 
transport dynamo. The period of the cycle in the theo- 
retical model is approximately inversely proportional to 
the amplitude of the meridional circulation [l(J EH • Pre- 
sumably the variations in the periods of past cycles were 
produced primarily by variations in the meridional cir- 
culation. Although we have direct measurements of the 
meridional circulation flow speed only during the last few 



years 



19l 2oj j . the periods of past cycles can be used 



to 



draw inferences about meridional circulation variations in 
the past 2lH24f. Such variations in meridional circula- 
tion also do cause variations in the strengths of different 
cycles in addition to variations in their periods. Sup- 
pose the meridional circulation has become weaker than 
usual. Then the period of the cycle will be longer and 
diffusion will have more time to act on the magnetic field, 
making the cycle weaker. If diffusion is assumed in the 



range 10 -10 cm s consistent with mixing length 
arguments [25 ] (which we do in our model), then this ef- 
fect overcomes the opposing effect of differential rotation 
also getting more time to generate more toroidal field 
and the cycles become weaker when the meridional cir- 
culation slows down 18|, |2lj, [24j . However, in the model 
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developed by the HAO group, the diffusion is taken to 
be about 50 times smaller than what we take [l^, Hfjjl, 
leading to the opposite effect of cycles getting stronger 
with slower meridional circulation due to the generation 
of more toroidal field by differential rotation over a longer 
time. Several arguments in favor of the higher diffusi vity 
used by us are summarized in Sect. 5 of Jiang et al. [33|. 

In order to model grand minima theoretically, we have 
to run our theoretical flux transport dynamo model with 
fluctuations in poloidal field generation and fluctuations 
in meridional circulation. Some studies have shown that 
large fluctuations in poloidal field generation mechanism 
or large fluctuations in the meridional circulation can 
force the dynamo into intermittencies resembling grand 
minima 27 -32|. Karak 21 1 found that the dynamo is 
driven into a grand minimum if the poloidal field at the 
end of a cycle and the meridional circulation at that time 
fall to sufficiently low values. Defining a grand minimum 
as absence of sunspots for at least 20 years (the same 
definition as used by Usoskin et al. Q in estimating the 
number of grand minima in the past 11,000 years), the 
shaded region of Fig. 1 indicates the combined values of 
poloidal field and meridional circulation at the end of a 
cycle necessary for forcing the dynamo into grand min- 
ima according to our theoretical model. The important 
question now is to estimate the probability of this hap- 
pening at the end of a cycle. We estimate this probability 
in the following way. 

Assuming that the inverse of the cycle period gives an 
approximate value of the meridional circulation ampli- 
tude during that cycle (as s ugg ested by theoretical flux 
transport dynamo models [id Hij|). Karak & Choud- 
huri [24J concluded that the meridional circulation has 
changed randomly in the past with a correlation time 
around 30-40 yr (also see the similar study based on the 
low order dynamo model 22J, |23[). Figure 2a shows a 



histogram of the estimated values of the meridional cir- 
culation amplitude during the last 28 cycles. The solid 
curve in Fig. 2a shows the Gaussian having the mean 
vq = 23 m s _1 and the standard deviation a v = 3.34 m 
s _1 calculated from the data presented in the histogram. 
We find that the Gaussian is a reasonable fit, although 
we do not have any data points lying way out in the 
Gaussian tail. Jiang et al. [33| (also see [15j]) pointed 
out on the basis of observational data that there is a 
good correlation between the poloidal field at the end of 
a cycle and the strength of the next cycle. Assuming 
a perfect correlation, we can use the strengths of fol- 
lowing cycles to obtain values 7 of poloidal field at the 
ends of previous cycles (scaled by taking their average 
value as the unit). Figure 2b is a histogram of the val- 
ues 7 of poloidal field obtained in this way. It should, 
however, be kept in mind that variations in the merid- 
ional circulation also contribute to fluctuations of cycle 
strengths [H, [2l|, 0, HH and fluctuations of poloidal field 
at the end of cycles obtained by our method (as shown 
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FIG. 1: The solid line shows the values of the meridional 
circulation amplitude vo and the poloidal field scale factor 7 
which produce grand minima of duration ~ 20 yrs. The pa- 
rameters lying in the shaded region produce grand minima 
of longer duration. The dashed curves are the contours of 
the joint probability P(j, «o), with the values of P(j, vo) (ex- 
cluding the constant pre-factor) given in the plot. The scale 
factor 7 is defined as the amplitude of poloidal field at the end 
of a cycle divided by its value in the absence of fluctuations. 
This figure is produced by using our theoretical solar dynamo 
model (details given in Karak [211]') . 
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FIG. 2: (a) Histogram of the meridional circulation amplitude 
derived from the observed periods of last 28 solar cycles. The 
solid curve is the Gaussian function with mean = 23 m s _1 
and standard deviation a v = 3.34 m s _1 . (b) Histogram of 
poloidal field scale factor derived from the peak sunspot num- 
ber of last 28 solar cycles. The solid curve is the Gaussian 
function with mean = 1 and standard deviation a 1 = 0.35. 
The probabilities plotted along the vertical axis in both (a) 
and (b) are obtained by dividing the number of data points 
in a data bin by the total number of data points and the 
horizontal width of the data bin. 



in Fig. 2b) is probably an over-estimate. The solid curve 
in Fig. 2b is the Gaussian having the standard deviation 
07 = 0.35 obtained from the data in the histogram (the 
mean is 1 by definition). We expect on general grounds 
the distribution of 7 to follow a Gaussian, although we 
find that the Gaussian fit in Fig. 2b is not as good a fit as 
in Fig. 2a. This is not surprising given that our data set 
comprising of only 28 cycles is quite small. If we assume 
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that the fluctuations in the amplitude vo of meridional 
circulation at the solar surface and the fluctuations in 7 
both follow Gaussian distributions, then we can draw one 
obvious inference. The joint probability that the poloidal 
field at the end of a cycle lies in the range 7, 7 + and 
the amplitude of meridional circulation at the same time 
lies in the range vq, vq + dvg is given by 



P('y,vo)d'jdv 



1 



a v y2n 



exp 
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dj dvQ 



In Fig. 1 we show various contours corresponding to 
different value of P(7, vq) (excluding the constant pre- 
factor). The probability that 7 and vq at the end of a 
cycle jointly lie within a certain area in Fig. 1 is easily 
obtained by integrating J P(y, vq) dj dvo over that area. 
On carrying out this integration over the shaded region 
in Fig. 1 that corresponds to conditions for producing 
grand minima, we find the value of the integral to be 
0.013 (i.e., 1.3%), remarkably close to the probability of 
occurrence of grand minima on the basis of observational 
data Table 1 gives the value of this integral for a 
few combinations of a v and 07. It may be noted that 
there would be contributions to the probability from the 
regions to the left and to the bottom of the shaded region 
in Fig. 1. We have checked that these contributions are 
negligible. 

Table 1. This table gives the values of the probability (in %) 
of the initiation of grand minima (given by J P(~y, vo) d-y dvo 
over the shaded region of Fig. 1) for different combinations of 
the standard deviations a v and 07. 





07 

0.25 0.35 0.46 


3.20 
a v 3.34 
3.50 


0.8 1.3 1.7 
0.9 1.3 1.9 
1.1 1.7 2.2 



To check whether grand minima really do occur in ac- 
cordance with the above simple probability estimate, we 
have carried out extensive simulations on the basis of our 
dynamo code 11, 35|. Karak [2lj changed the values of 
a few parameters and our present simulations are based 
on this model. We introduce fluctuations in poloidal field 
gen eration by the method proposed by Choudhuri et al. 
17| . At the end of every cycle, we multiply the poloidal 
field above 0.8i?© by the random number 7 obeying the 
Gaussian distribution shown in Fig. 2b. This procedure 
introduces fluctuations in the poloidal field generated in 
the last cycle lying in the upper portions of the convec- 
tion zone, whereas any poloidal field produced at the 
earlier cycles lying at the bottom of the convection zone 



remains unchanged. While this procedure introduces a 
momentary discontinuity in the field lines at depth O.8-R0 
(see Fig. 1 in Jiang et al. [ill), this discontinuity disap- 
pears soon and does not cause any problem. To introduce 
fluctuations in meridional circulation, we change its am- 
plitude randomly after every 30 yr such that the ampli- 
tudes obey the Gaussian distribution shown in Fig. 2a. 
On running the code for 11,000 yr in this way, we ob- 
tain the number of grand minima in the range 24-30 for 
different realizations of randomly generated meridional 
circulation and 7 — remarkably close to the observational 
finding of 27 grand minima in the last 11,000 yr. Another 
important result is that we find the Sun to spend about 
10-15% of the time in a grand minimum state, which is 
very close to 17% found in the observational study [3]. 
Fig. 3 shows the durations of grand minima and their 
times of occurrence for a particular run. We ourselves 
have been amazed that the observational data are repro- 
duced so well on running our code with the simple incor- 
poration of the fluctuations suggested by the histograms 
in Fig. 2, without having to change any parameters of the 
dynamo model compared to our previous work. Fig. 4 is 
a sample plot showing how the sunspot number varied 
for a typical 1000 years during which two grand minima 
occurred. Histograms of the duration of grand minima 
and waiting time between grand minima are shown in 
Fig. 5. These histograms are to be compared with Figs. 7 
and 6 in Usoskin et al. In the simulations presented 
above, we have changed meridional circulation abruptly 
after every 30 yr. To check how the above results change 
with this coherence time r, we have done several simula- 
tions with different r. We see that even when r ~ 15 yr, 
we get around 15-20 grand minima. However, when r 
is less than this value, the number of grand minima gets 
very much reduced. Another important point: instead of 
changing meridional circulation abruptly at a time, if we 
change it smoothly in few years, then also the results do 
not change significantly. 

When sunspots are absent, the Babcock-Leighton pro- 
cess for the generation of poloidal field cannot take place. 
Presumably, during a grand minimum, the poloidal field 
has to be generated by the a-effect originally proposed 
by Parker j3(| and Steenbeck, Krause & Radler [37j |. 
It is possible that an a-effect coexists along with the 
Babcock-Leighton mechanism all the time, although its 
nature and spatial distribution (even its sign) is com- 
pletely unclear at this time. In view of this uncertainty, 
our results presented in Figs. 3 and 4 are obtained by 
using the same form of a at all times. Once the Sun is 
pushed into a grand minimum, it comes out of the grand 
minimum in a time of the order of dynamo growth time, 
as discussed in detail by Choudhuri & Karak 29]. One 
intriguing fact to note is that we get durations of grand 
minima similar to their observed values when we keep 
using the same a throughout the grand minima. 

We conclude that the irregularities of solar cycles in- 
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FIG. 3: The durations of grand minima indicated by vertical bars at their times of occurrence in a 11,000 yr simulation. This 
is the result of a particular realization of random fluctuations that produced 28 grand minima. 
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FIG. 4: The plot of sunspot number against time for a typical 
1000 years. The two grand minima around 1100 and 1500 
years can be seen in Fig. 3 as well. 
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FIG. 5: Left panel shows the distribution of the durations of 
the grand minima and the right panel shows the distribution 
of the waiting times between the grand minima. 



eluding the grand minima are produced by fluctuations 
in poloidal field generation and in meridional circulation. 
Assuming these fluctuations to obey Gaussian distribu- 
tions, we obtain the basic parameters of the Gaussians 
from the distribution of the peak sunspot numbers and 
the durations of the last 28 solar cycles. On running our 
code with such fluctuations, we find that our theoretical 
dynamo model produces grand minima at a frequency 
remarkably close to what is found in the data for the last 
11,000 yr. 
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